Symmetries and Reversing Symmetries 
of Toral Automorphisms 



Michael Baake 1 ) and John A. G. Roberts 2) 

1) Institut fur Theoretische Physik, Universitat Tubingen, 
Auf der Morgenstelle 14, 72076 Tubingen, Germany 

2) Department of Mathematics, LaTrobe University, 
Bundoora, Victoria 3083, Australia 



Abstract 

Toral automorphisms, represented by unimodular integer matri- 
ces, are investigated with respect to their symmetries and reversing 
symmetries. We characterize the symmetry groups of GL(n, Z) matri- 
ces with simple spectrum through their connection with unit groups 
in orders of algebraic number fields. For the question of reversibility, 
we derive necessary conditions in terms of the characteristic polyno- 
mial and the polynomial invariants. We also briefly discuss extensions 
to (reversing) symmetries within affine transformations, to PGL(n, Z) 
matrices, and to the more general setting of integer matrices beyond 
the unimodular ones. 



Introduction 

Unimodular integer matrices induce interesting dynamical systems on the 
torus, such as Arnold's famous cat map 0, Ch. 1, Ex. 1.16]. This is an 



example of a hyperbolic dynamical system that is ergodic and mixing |27| , 
and also a topological Anosov system. Therefore, with a suitable metric, 
it makes the 2-torus into a Smale space, see Thm. 1.2.9] for details. 



Although induced from a linear system of ambient space, the dynamics on 
the torus is rather complicated, and these systems serve as model systems in 
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symbolic dynamics and in many applications. Very recently, also cat maps 
on the 4-torus (and their quantizations) have begun to be studied |TT| . 

Hyperbolic toral automorphisms also play a prominent role in the theory of 
quasicrystals through their appearance as inflation symmetries, see || |8j and 
references therein. In particular, in the one-dimensional case, these symme- 
tries give rise to interesting non-linear dynamical systems called trace maps 
J?], [32|, [L4| that have extensively been used to study the physical properties 
of one-dimensional quasicrystals. 

It is always helpful to know the symmetries and reversing symmetries of 
a dynamical system [[53|, |22| , and it was perhaps a little surprising that in 
the case of GL(2, Z) and PGL(2, Z) a rather complete classification could 
be given, see flO| for a detailed account or |9|, || for a summary. The main 
difficulty when the matrix entries are restricted to integers is that one can 
no longer refer to the usual normal forms of matrices over C or R, but has to 
use discrete methods instead. Fortunately, there is a strong connection with 
algebraic number theory, see [37| for an introduction, and this connection is 
certainly not restricted to the 2D situation. 

It is thus the aim of this article to extend the results of our earlier article 
TU| to the setting of matrices in GL(n, Z). The answers will be less complete 



and also less explicit, but the connection to unit groups in orders of algebraic 
number fields is still strong enough to give quite a number of useful and 
general results, both on symmetries and reversing symmetries. From a purely 
algebraic point of view, the results derived below are actually rather straight- 
forward. However, these results, and the methods used to derive them, are 
not at all common in the dynamical systems community. Therefore, this 
article is also intended to introduce some of these techniques, and we try to 
spell out the details or give rather precise references at least. Furthermore, 
as with our article [ITJ], the results of this paper have relevance to both the 
dynamics community (e.g. the dynamics of hyperbolic toral automorphisms 
generated by (symplectic) SL(4, Z) matrices [01J]) and to the quasicrystal 
community (e.g. inflation symmetries of planar point sets projected from 4D 
lattices, where the symmetries are generated by GL(4, Z) matrices || |§). 

The article is organized as follows. We start with a section on the back- 
ground material, including the group theoretic setup we use and a recollection 
of those results from algebraic number theory that we will need later on. Sec- 
tion is the main part of this article. Here, we derive the structure of the 
symmetry group of toral automorphisms with simple spectrum and discuss 
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reversibility. Section ^ extends the set of possible symmetries to affine trans- 
formations and summarizes the analogous problem for projective matrices. 
It also discusses the extension of (reversing) symmetries to matrices that are 
no longer unimodular. 



1 Setting the scene 

In this Section, we explain in more detail what we mean by symmetries 
and reversing symmetries, and we also recall some results from algebra and 
algebraic number theory that we will need. 



1.1 Symmetries and reversing symmetries 

For a general setting, consider some (topological) space Q, and let Aut(O) be 
its group of homeomorphisms or, more generally, a subgroup of homeomor- 
phisms of Q which preserve some additional structure of Q. Consider now an 
element F 6 Aut(f2) which, by definition, is invertible. Then, the group 

S(F) := {G e Aut(fi) | GoF = FoG} (1) 

is called the symmetry group of F in Aut(fi). In group theory, it is called the 
centralizer of F in Aut(O), denoted by cent Aut ^ (F) . This group certainly 
contains all powers of F, but often more. 

Quite frequently, one is also interested in mappings R G Aut(f2) that 
conjugate F into its inverse, 

RoFoR- 1 = F~K (2) 

Such R is called a reversing symmetry of F, and when such an R exists, we 
call F reversible. We will, in general, not use different symbols for symmetries 
and reversing symmetries from now on, because together they form a group, 

U{F) := {G e Aut(fi) \ GoFoG- 1 = F^ 1 } , (3) 



the so-called reversing symmetry group of F, see |2(J for details. If (F) 
denotes the group generated by F, TZ(F) is a subgroup of the normalizer of 
(F) in Aut(fi). 

There are two possibilities: either 1Z(F) = S(F) (if F is an involution or 
if it has no reversing symmetry) or TZ(F) is a C2-extension (the cyclic group 
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of order 2) of S(F) which means that S(F) is a normal subgroup of TZ(F) 
and the factor group is 

K{F)/S{F) ~ C 2 . (4) 
The underlying algebraic structure has fairly strong consequences. One is 



that reversing symmetries cannot be of odd order pD |, another one is the 
following product structure [|It], Lemma 2]. 

Fact 1 If F (with F 2 ^ Id) has an involutory reversing symmetry R, the 
reversing symmetry group of F is given by 

H(F) = S(F) x s C 2 , (5) 

i.e. it is a semi-direct product^ □ 

We can say more about the structure of 7Z(F) if we restrict the possibilities 
for S(F), e.g. if we assume that S(F) ~ or S(F) ~ x C 2 with the C 2 
being a subgroup of the centre of Aut(fi), compare also [|16[ for some group 



theoretic discussion. This situation will appear frequently below. 

1.2 (Reversing) symmetries of powers of a mapping 

In what follows, we summarize some of the concepts and results of Ref. ET 



and, in particular, Ref. pOfl . It may happen that some power of F has more 
symmetries than F itself (we shall see examples later on), i.e. S(F k ) (for some 
k > 1) is larger than S(F) which is contained as a subgroup. The analogous 
possibility exists for TZ(F k ) versus H(F). If such a situation occurs, we say 
that F possesses additional (reversing) /c-symmetries. Let us make this a 
little more precise. 

It is trivial that mappings F of finite order (with F k = Id, say) possess 
the entire group Aut(fi) as ^-symmetry group. Let us thus concentrate on 
mappings F G Aut(fi) of infinite order. We denote by Sao(F) the set of 
automorphisms that commute with some positive power of F . This set can 
be seen as the inductive limit of S(F k ) as k — > oo, with divisibility as partial 
order on N, and Soo(F) is thus a subgroup of Aut(f2). Let # F (G) denote the 
minimal k such that G o F k = F k o G. Then 

Soo(F) = {Ge Aut(fi) | # F (G) < oo } . (6) 



1 We use N x s H for the semi-direct product of two groups N and H, with N being the 
normal subgroup. 
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Of course it may happen that #f{G) e 1 on Soo(F) which means that 
no power of F has additional symmetries. On the other hand, #f(G) might 
be larger than one in which case we call G a genuine or true fc-symmetry. 
G is a truef] /c-symmetry of F if and only if the mapping G i— > F o G o F" 1 
generates a proper fc-cycle. We shall meet this phenomenon later on. 

Quite similarly, one defines reversing fc-symmetries and their orbit struc- 
ture [p0|| , but we will not expand on that here. 



1.3 Some recollections from algebraic number theory 

Much of what we state and prove below can be seen as an application of 
several well-known results from algebraic number theory. The starting point 
is the connection between algebraic number theory and integral matrices, see 



37J for an introduction. 

To fix notation, let Mat(n, Z) denote the ring of integer^ n x n-matrices. 
An element M of it is called unimodular if det(M) = ±1, and the subset 
of all unimodular matrices forms the group GL(n, Z). For the characteristic 
■polynomial of a matrix M, we will use the convention 

n 

P{x) := det(xl-M) = JJ(ar-Ai) (7) 

8=1 

where Ai,...,A n denote the eigenvalues of M. With this convention, P(x) 
is monic, i.e. its leading coefficient is 1. If M is an integer matrix, P(x) 
has integer coefficients only, so all eigenvalues of M are algebraic integers. 
Conversely, the set of algebraic integers, which we denote by A, consists of 
all numbers that appear as roots of monic integer polynomials. 

To show the intimate relation more clearly, let us recall the following 
property (see item (b) on p. 306 of 



Fact 2 Let P(x) = x n + a n _ x x n ~ x + . . . + a x x + a Q be a monic polynomial 
with integer coefficients a i that is irreducible over Z. Let a be any root of 
it, and A an integer matrix that has P(x) as its characteristic polynomial. 
Then, the rings Z[ot] and 7h\A\ are isomorphic. □ 



2 Although the distinction between true and other fc-symmetries is necessary in general, 
we shall usually drop the attribute "true" whenever misunderstandings are unlikely. 

3 Here, and in what follows, integer means rational integer, i.e. an integer in Q. Other 
kinds of integers, such as algebraic, will be specified explicitly. 
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We write Z[x] for the ring of polynomials in x with coefficients in Z, see 
S3] , p. 90] for details. Clearly, Fact [2] also extends to the isomorphism of the 



rings Q[a] and Q[A]. For background material on polynomial rings, we refer 



to [23 , Ch. IV]. Let us only add that a polynomial in Z[x] is irreducible over 
Z if and only if it is irreducible over Q, see [p3j Thm. IV. 2. 3]. 

Let P(x) be any (i.e. not necessarily irreducible) monic integer polynomial 
of degree n. In general, there are many different matrices A which have P(x) 
as their characteristic polynomial, and even different matrix classes (we say 
that A,B G Mat(n, Z) belong to the same matrix class if they are conjugate 
by a GL(n, Z) matrix C, i.e. A = CBC~ l ). Let us recall the following helpful 
result on the number of matrix classes, see [[57], Thm. 5] for details: 

Fact 3 Let P(x) be a monic polynomial of order n with integer coefficients 
that is irreducible over 7L, and let a be any of its roots. Then the number of 
matrix classes generated by matrices A G Mat(n, Z) with P{A) = equals 
the number of ideal classes (or class number, for short) of the order Z[a]. 
In particular, this class number is finite, and it is larger than or equal to the 
class number of the maximal order O max o/Q(a). □ 

Let us explain some of the terms used here. If a is an algebraic number, 
Q(a) denotes^] the smallest field extension of the rationals that contains a. Its 
degree, n, is the degree of the irreducible monic integer polynomial that has 
a as its root. The set O max := Q(a) fl A is the ring of (algebraic) integers 
in Q(a), and is called its maximal order. More generally, a subring O of 
C^max is called an order, if it contains 1 and if its rational span, QO, is all of 
Q(a). A subset of O is called an ideal if it is both a Z-module (i.e. closed 
under addition and subtraction) and closed under multiplication by arbitrary 
numbers from O. The ideals come in classes that are naturally connected to 



the matrix classes introduced above, see O, |37j for further details. 



An element e G Q(a) is called a unit (or, more precisely^, a unit in C max ) 
if both e and its inverse, e -1 , are algebraic integers and hence are in max . 
This happens if and only if the corresponding matrix is unimodular, i.e. if any 
monic integer polynomial P(x) that has e as a root has coefficient a Q = ±1. 
So, matrices in GL(n, Z) and units in algebraic number fields are two facets 
of the same coin. The units of O max form a group under multiplication, 



4 Note the difference between the meaning of Q[a] and 

5 This distinction is useful if orders other than O max appear. 
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denoted by C^ ax in the sequel. Similarly, if O C O m3iX is any order of Q(a), 
we write O x for its group of units, which is then a subgroup of C^ax- 

Independently of whether the monic polynomial P(x) is irreducible or not, 
there is always at least one matrix which has characteristic polynomial P(x), 
namely the so-called (left) companion matrix: 



( o 





A (i) 



1 





V 





-On 



J 



Consequently, the class number is > 1 (if P(x) is irreducible, the companion 



matrix actually corresponds to the principal ideal class, see [p7| Thm. 9]). 
Another obvious choice, always belonging to the same matrix class, is the 
(right) companion matrix, A^ r \ obtained from by reflection in both 
diagonal and anti-diagonal, i.e. 



A (r) = RA W R 



with the involution 



R 



(° 

V 



(9) 



(10) 



1 °J 

This matrix will reappear several times in what follows. 



2 GL(n, Z) matrices and toral automorphisms 

Let us generally assume that n > 2. The toral automorphisms of the n-torus 
T n := M n /Z n can be represented by the unimodular n x n-matrices with 
integer coefficients which form the group GL(n, Z). It now plays the role of 



Aut(fi) from Section Note that the elements of GL(n, Z) preserve the 



linear structure of the torus. 



2.1 Symmetries 



The first thing we will look at, given a toral automorphism M G GL(n, Z), 
is its symmetry group within the class of toral automorphisms. So, we want 
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to determine the centralizer of M in GL(n, Z), 

S(M) = cent GL(niZ) (M) = {G G GL(n, Z) | MG = GM } . (11) 

To be more precise, we are mainly interested in the structure of the symmetry 
group rather than in explicit sets of generators and relations. This is invariant 
under conjugation, i.e. if we know it for an element M, we also know it for 
any other element of the form BMB^ 1 because 

S(BMB^ 1 ) = BS{M)B~ l . (12) 

A given integer matrix M G GL(n, Z) determines its characteristic polyno- 
mial (0) which is monic and has integer coefficients, so its roots are algebraic 
integers. Now, two principal situations can occur for the characteristic poly- 
nomial: it is either reducible over Z (which happens if and only if at least 
one eigenvalue of M is an algebraic integer of degree less than n) or it is ir- 
reducible. In the latter case, since we are working over the field Q, we know 
that the roots must be pairwise distinct. So we have 

Fact 4 Let M be an integer matrix with irreducible characteristic polyno- 
mial. Then M is simple and hence diagonalizable over C. □ 



Here, M is called simple if it has no repeated eigenvalues (which is also called 
separable elsewhere). This case will be dealt with completely. 

If the characteristic polynomial is reducible, the matrix can still be simple, 
and we will see the general answer for this case, too. Beyond that, M can 
either be semi-simple (i.e. diagonalizable over C) or not, and we will not 
be able to say much about this case. This is really not surprising, as this 
situation is closely related to the rather difficult classification problem of 
crystallographic point groups, see [1T2] for answers in dimensions < 4 and 
p9| for a recent survey. 

Let us now state one further prerequisite for tackling the symmetry ques- 
tion. In view of later extensions, we do this in slightly more generality 
than needed in the present Section. Recall that an n x n-matrix M, acting 
on a vector space V, is called cyclic, if a vector v G V exists such that 
{v, Mv, M 2 v, ... , M n ~ x v } is a basis of V. Also, the monic polynomial Q of 
minimal degree that annihilates M, i.e. Q(M) = 0, is called the minimal 
polynomial of M. By the Cayley-Hamilton Theorem, it always is a factor of 
the characteristic polynomial of M. 
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Fact 5 Let M e Mat(n, Q) 6e a rational matrix, with characteristic polyno- 
mial P(x) and minimal polynomial Q(x). Then the following assertions are 
equivalent. 

(a) The matrix M is cyclic. 

(b) The degree of Q(x) is n. 

(c) P{x) = Q{x). 

(d) G G Mat(n, Q) commutes with M <==^ G 6 Q[M]. 

PROOF: A convenient source is [18, Ch. III]. The equivalence of statements 



(a) - (c) is a consequence of Thm. III. 2. The equivalence of (a) with (d) 
follows from Thm. III. 17, and the Corollary following it, together with the 
Corollary of Ch. III. 17. Alternatively, see |, Cor. 5.5.16]. □ 

Lemma 1 Let M e GL(n, Z) have a characteristic polynomial P(x) that is 
irreducible over Z, and let X be a root of P(x). Then the centralizer of M in 
GL(n, Z) is isomorphic to a subgroup of finite index of the unit group in the 
ring of integers C max of the algebraic number field Q(A) . 

Proof: By assumption, P(x) is also the minimal polynomial of M, so any 
GL(n, Z)-matrix which commutes with M is, by Fact a polynomial in M 
with rational coefficients. Consequently, S(M) is isomorphic to a subset of 
Q[M] that forms a group under (matrix) multiplication. So, we have to 
analyze Q[M] to find out what this group is. 

Let (P(x)) denote the ideal in Q[x] generated by our polynomial P(x). 



Then Q[x}/(P(x)) ~ Q[A], see p|, p. 224], and Q[A] ~ Q[M], by Fact|resp. 
the remark following it. Since A is algebraic over Q and P(x) is irreducible, 
we know by [^, Prop. V.1.4] that Q[A] = Q(A) is an algebraic number field, 
of degree n over Q. Under the isomorphism Q[M] ~ Q(A), GL(n, Z)-matrices 
correspond to units in Q(A), hence S(M) must be isomorphic to a subgroup 
of Cmax; the unit group of the maximal order of Q(A). 

Observe that every matrix in Z[M] commutes with M, in particular those 
of Z[M] n GL(n,Z), which form a subgroup of S(M). But Z[M] ~ Z[A] 
means that this subgroup is isomorphic to the unit group Z[A] X . So, if we 
identify S(M) with its image in Q(A) under the isomorphism, it is sandwiched 
between Z[A] X and 0^ ax . Note that Z[A] C O max is an order. 
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Finally, recall that the unit group of an order O is a finitely generated 
Abelian group, and that it is always of maximal rank, see [T2^, Thm. 1.12.12] 
or [pi Sec. 4, Thm. 5], i.e. its rank equals that of the unit group of the 
maximal order O max . In particular, the group-subgroup index [O max : Z,[A] X ] 
is finite, and S(M) must then also be of finite index in max . □ 

Let us comment on this result. First of all, it does not matter which root 
Aj of P(x) we choose, as all the n (possibly different) realizations Q(Aj) are 
mutually isomorphic, and so are their unit groups. Explicit isomorphisms 
are given by the elements of the Galois group of the splitting field /C = 
Q(A 1 , ... , A n ) of P(x), see [|^, Ch. VI. 2] for details. Note also that, in general, 
Q(A) will be a true subfield of the splitting field K - so, it is really the unit 
group of Q(A) that matters, and not the unit group of /C. 

Another way to view the result, in a more matrix oriented way (and similar 
to our approach in |K|), is to look at the diagonalization of M, 

UMU- 1 = diag(A 1 ,...,AJ. (13) 

Here, £7 _1 can be arranged to have its j-th column in the field Q(A ■) because 
one can solve the corresponding eigenvector equation in the smallest field 
extension of Q that contains A,-. In fact, we only have to do this for the 
first column — the others are then obtained by applying appropriate Galois 
automorphisms to the first one. 

Any other matrix G G GL(n, Z) with [G, M] = GM - MG = must now 
also fulfil 

[UGU~\UMU- 1 ] = UiG^Mp- 1 = 0. 

But only diagonal matrices can commute with diag(A 1; ... , A n ) = UMU^ 1 
because the eigenvalues are pairwise distinct. So, we must have 

UGU' 1 = diag(/i 1; ... ,//J , 

with all fii G Q(Aj) units. They are, however, not independent but obtained 
from one another by the same set of Galois automorphisms that were used 
to link the columns of the matrix U~ l , which is why we get the result. 

Lemma [3] raises the question: What is the unit group of the maximal 
order in K = Q(A) ? The answer is given by Dirichlet's unit theorem, see 
[|T3| , Sec. 11. C] or p. 334]. Group the roots of the irreducible polynomial 
P(x) into n 1 real roots and n 2 pairs of complex conjugate roots, so that 
n = n l + 2n 2 . (In other words: we have n l real and n 2 pairs of complex 
conjugate realizations of the abstract number field 
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Fact 6 Let X be an algebraic number of degree n = n l + 2n 2 , with n 1 and 
n 2 as described above. Then, the units in the maximal order O max of the 
algebraic number field K = Q(A) form the group 

E(K) = max - Tx Z"!^- 1 (14) 

where T = (9 max H {roots of unity} is a finite Abelian group and cyclic. □ 

In particular, this means that T, which is also called the torsion subgroup 
of E(K), is generated by one element. In many cases below, we will simply 
find T ~ C*2. Combining now Lemma [I] with Fact we immediately obtain 

Proposition 1 Under the assumptions of Lemma [I], the symmetry group 
S(M) C GL(n,Z) is a subgroup of E(K) of (|H|) of maximal rank, i.e. we 
have 

<S(M) ~ T' x Z"!^" 1 
where T' is a subgroup of the torsion group T as it appears in ( |!H| ) . □ 

Note that Proposition [TJ does not imply that the torsion-free parts of S(M) 
and E(K) are the same, only that they are isomorphic. In fact, a typical 
situation will be that they are different in the sense that E(K) is generated 
by the fundamental units, but S(M) only by suitable powers thereof. 

What, in turn, can we say about the torsion group T' in Proposition |l|? 
Whenever the characteristic polynomial P(x) of M 6 GL(n, Z) is irreducible 
and has at least one real root (e.g. if n is odd), a say, then K = Q(ct) is 
real, and Q(a) H S 1 = {±1}, where S 1 is the unit circle. Consequently, the 
torsion subgroup of E(K) in this case is T = {±1} ~ C 2 . Since a toral 
automorphism always commutes with ±1, we obtain 

Corollary 1 //, under the assumptions of Lemma |l|, one root of the irre- 
ducible polynomial P(x) is real, the torsion group in Proposition [I] is T' ~ 
Ln particular, this is the case whenever the degree of P(x) is odd. □ 

Let us look at two examples in GL(3,Z), namely 

/o o i\ n 

M 1 = 1 and M 2 = 1 1 I , (15) 
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which are taken from |2J], Eqs. (5.21) and (5.3)]. They have been studied 
thoroughly in the context of inflation generated one- dimensional quasicrys- 
tals with a cubic irrationality as inflation factor. We have det(M 1 ) = 1, 
det(M 2 ) = —1, and the characteristic polynomials are Pi(x) = x 3 — x 2 — 1 
and P2{x) = x 3 — 2x 2 — x + 1, both irreducible over Z. Both matrices are 
hyperbolic, and the largest eigenvalue in each case is a Pisot-Vijayaraghavan 
number, i.e. an algebraic integer > 1 all algebraic conjugates of which lie 
inside the unit circle, see [55 for details. 



Now, Mi has one real and a pair of complex conjugate roots, so our above 
results lead to S(M{) ~ C2 x Z, where the infinite cyclic group is actually 
generated by M\ itself because its real root is a fundamental unit. M2, in 
turn, has three real roots, and we thus get S(M2) — C2 x Z 2 . As generators 
of Z 2 , one may choose M 2 and M' 2 where 



M' 2 



which can be checked explicitly. This example is of relevance in connection 
with planar quasicrystals with sevenfold symmetry, see [21, Sec. 5.2] for de- 




tails, where the three-dimensional toral automorphism M 2 shows up in the 
cut and project description of special directions in the quasicrystal. Other 
examples related to planar quasicrystals with 8-, 10- and 12-fold symmetry, 
in which the torsion subgroup T' of Proposition [l] is different from C2, will 



be given in Section 2.2 



Let us now return to the general discussion and extend the previous results 
to the case that M is simple, and hence diagonalizable (over C) with pairwise 
different eigenvalues. Since diagonal matrices with pairwise different entries 
only commute with diagonal matrices, we have: 

Corollary 2 If M G GL(n, Z) is simple, <S(M) C GL(n,Z) is Abelian. □ 

Note that the converse is not true: even if M is only semi-simple, or not even 
that (i.e. not diagonalizable), S(M) can still be Abelian, e.g. if M observes 
the conditions of Fact E. As far as we are aware, not even the Abelian 



subgroups of GL(n, Z) are fully classified, see |26|, |2j| and references given 
there for background material. 

Let P{x) be the characteristic polynomial of a simple matrix M. If it is 
reducible over Z, it factorizes as P{x) = rii=i-^( x ) irreducible monic 
polynomials P{(x). 
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Theorem 1 Let M G GL(n, Z) 6e simple and let its characteristic polyno- 
mial be P(x) = Y[i=i Pi{ x )j with Pi{x) irreducible over Z. Then, the symme- 
try group of M, S(M) C GL(n, Z) ; is a finitely generated Abelian group of 
the form 

S(M) = T x U (16) 

where T is a finite Abelian group of even order, with at most £ generators. 
Furthermore, if the irreducible component Pi(x) has njj 1 real roots and 
pairs of complex conjugate roots, the rank r of the free Abelian group in ( JTB| ) 
is given by 

t 

r = Y^{nf +nf -1). (17) 

i=i 

PROOF: Since M is simple, the degree of its minimal polynomial is n and 
Fact [5] tells us that cent Mat ^ n ^ = Q[M}. As P{x) has no repeated factors 
(so that Q[M] contains no radicals), we get, by []T8| , Thm. III. 4], 

Q[M] ~ Q[aJ © ... © Q[a e ] (18) 

where a i is any root^ of Pi(x), for 1 < i < t. 

Under the assumptions made, each Q[aJ = Q(aJ is a field. Since a 
GL(n, Z)-matrix in Q[M] will correspond to a unit in each of the Q(ojj), we 
can now apply Lemma [I] to each component, giving ([H]) as the direct sum 
of i unit groups. The rank in ([17|) follows now from Proposition [l]. 

The torsion part T is a finitely generated Abelian group, with (at most) 
one generator per irreducible component of P(x), of which there are I. 
Clearly, S(M) always contains the elements ±1, so {±1} — Ci is a sub- 
group of T. The order of T is then divisible by 2, hence even. □ 



Although Theorem does not give the general answer to the question for 
the symmetry group S(M), it certainly gives the generic answer, because 
the property of M having simple spectrum is generic. But what about the 
remaining cases? Without further elaborating on this, let us summarize a 
few aspects and otherwise refer to the literature ^6, [29| for a summary of 
methods to actually determine the precise centralizer. 

6 Note that the oi i have pairwise different minimal polynomials by assumption, but that 
Q[aJ ~ Q[ctj] for i ^ j is still possible. 
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If M G GL(n, Z) is semi-simple, but not simple, its characteristic polyno- 
mial contains a square, and whether or not S(M) is still Abelian (and then 
of the above form) depends on whether or not the minimal polynomial of M 
has degree n, see Fact [5]. Note, in particular, that the following situation can 
emerge. If P(x) has a repeated factor, but the corresponding matrix M is 
a block matrix, then the two blocks giving the same factor of P(x) can still 
be inequivalent, if the corresponding class number is larger than one which 
equals the number of different matrix classes, see Fact |3|. 

If M is not even semi-simple, things get even more involved. We can still 
have Abelian symmetry groups, e.g. if M is a Jordan block such as (J }), 
compare the results of [ID, Sec. 2.1.2] on parabolic automorphisms of T 2 . 
Clearly, this also follows from Fact [5]: if M G GL(n, Z) is conjugate to a 
single Jordan block, its minimal polynomial has degree n and all GL(n, Z)- 
matrices which commute with M are in Q[M]. This remains true if such a 
block occurs in a matrix that otherwise has simple spectrum disjoint from 1. 



Corollary 3 Let M G GL(n, Z). If the minimal polynomial of M has degree 
n, then S{M) C GL(n, Z) is Abelian. □ 

In a wider setting for symmetries, a stronger statement can be formulated, 
see Proposition [5] below and the comments following it. 

The general classification, however, and the non-Abelian cases in particu- 
lar, gets increasingly difficult with growing n and has been completed only for 
small n, see [^6| and references given there. Nevertheless, for any given M, 



the centralizer can be determined explicitly by means of various algorithmic 
program packages. 

Let us, at the end of this part and before we illustrate some of the above 
results by further examples, give a particular case of one matrix written as 
a polynomial of another. 

Fact 7 Let K be a field and M G GL(n, K) be an invertible matrix with 
characteristic polynomial P(x) = ^"=o a ^ x£ > where a n — 1 and a ^ 0. 
Then, the inverse matrix is given by 

n—l 

M- 1 = --$>+i M£ - 
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PROOF: Observe that P(M) = from the Cayley-Hamilton Theorem. The 
verification of M~ l M = 1 is then a straight-forward calculation. □ 

Before we discuss the connection of our approach to quasicrystallography 
in a separate Section, let us illustrate Theorem [I] with a recent example of a 
4D cat map taken from Eq. 3.21], namely 



M 



/o 



1 

V° 1 



1 




2 
1 



o\ 

-1 
1 

2 / 



(19) 



The characteristic polynomial is P(x) 
reducible over Z and splits as 



x 



4x 3 + 5x 2 — Ax + 1 which is 



P{x) = P 1 (x)P 2 (x) = {x 2 -3x + l)(x 2 -x+ 1) 

into Z-irreducible polynomials. Since P\ has two real and P 2 one pair of 
complex conjugate roots, Theorem [l| gives S(M) ~ T'xZ, with T' a subgroup 
of T = C2 x Cg. Also, since no root of Pi is a fundamental unit of the 
corresponding maximal order (which is Z[r] with r = (1 + y / 5")/2), the 
generator of the infinite cyclic group in S(M) could still differ from M. 

To determine the details, one easily checks that the most general matrix 
to commute with M is 

-c -d \ 

—d —c 
a + 2c + d a + c + 2d 
a + c + 2d a + 2c + dj 

A necessary condition for G to be in GL(4, Z) is then a, b,c,d G Z. This 
allows to exclude the existence of a root of M in S(M), and also no element 
of third order is possible. So we obtain 

S(M) = C 2 xC 2 x (M) . 

We will revisit this example below in the context of reversibility. 



G 



(a 
b 

c 

yd 



b 

a 
d 

c 



2.2 Three examples from planar quasicrystallography 

Planar tilings with 8-, 10- and 12-fold symmetry play an important role in 
the description of so-called quasicrystalline T-phases, see for background 
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material. They are of interest also in the present context because hyperbolic 
toral automorphisms show up through their inflation symmetry. 
For the 8-fold case, consider the polynomial 



P(x) 



x 4 + l 



(20) 



which has £, £ 3 , £ 5 , and £ 7 as roots, £ = e 2?r4//8 , which are primitive. So, P(x) 
is irreducible over Z, and Lemma [I] and Proposition [I] apply. In fact, Q(£) 



here is a cyclotomic field [38 with class number one, maximal order Z[£] and 
unit group Z[(] x ~C 8 xZ. 

If we denote the actual matrices that represent the generators for the 
groups C 8 and Z by M and G, respectively, it is natural to take the companion 
matrix of P(x) for M and to choose G accordingly, resulting in 



M 



(o 


1 





o\ 








1 














1 


[- 1 








V 



G 



( 


1 


l 





-A 




1 


l 


i 










l 


i 


l 


V 


-l 





i 





(21) 



By construction, M is a matrix of order 8. So, S(M) = (M,G) ~ C 8 x Z. 
What is more, anticipating the next Section, M turns out to be reversible 
with the matrix R of (|H]) as reversing symmetry. Using Fact [I] and observing 
[G, R] = 0, this means that TZ(M) = (M, G, R) ~ D 8 x Z. This, together 
with two similar examples, is summarized in Table [l|. 

Note that the case of 12-fold symmetry is more complicated because the 
fundamental unit in Z[£], for £ = e 2?ri//12 , is the square root of (2 + 
and hence not a simple homothety. This means that the representing matrix 
does not commute with R. The reversing symmetry group of this case, 
[Cyi x Z) x s C2, does contain a subgroup of the form D\ 2 x Z though - 
it is generated by M, G' = M~ 1 G 2 and R, where G' corresponds to the 
non-fundamental unit 2 + y/3. 



2.3 Reversibility 



The examples of Section |2^ were reversible, i.e. they fulfilled GMG~ l = M~ l 
for some G G GL(4, Z), in particular for the involution R e GL(4, Z) of fllCf). 
It is easy to check that this is also true for M of flHf). However, as we will 



see below, neither of the examples of (|i5|) are reversible in GL(3, Z 
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P{x) 

M 
(M) 

G 

S{M) 
K{M) 



( 1 
1 


V- 



8-fold 
x 4 + 1 

/ 1 o o\ 

10 
1 
\-\ 0/ 

C 8 

i o -A 

1 1 

1 1 1 

10 1 1 j 

C 8 xZ 

D R xZ 



10-fold 



x A + X 3 + x 2 + X 



1 



/0 



V- 

/ i 

-i 



1 



1 








1 





-1 -1 

c 5 

1 

-1 
-1 

1 

Cio x Z 
D 10 x Z 



+ 1 

o\ 



1 

-v 

i\ 



-1 



12-fold 

x 2 + 1 

1 o\ 

1 

1 
\-l 1 0/ 

Cl2 

1 o o\ 

1 1 
1 1 
1 1 1/ 

Cl2 x Z 

(Ci 2 x Z) x s C 2 



/o 






/ 1 






Table 1: Symmetries and reversing symmetries for three examples from qua- 
sicrystallography. The symmetry group is S(M) = (±M,G), and, similarly, 
K(M) = (±M, G, R), with R as in Eq. (|10|) for all three examples. 



In this Section, we are concerned with determining when reversibility can 
occur in GL(n, Z), and what we can say about the nature of the reversing 
symmetry G G GL(n, Z), e.g. whether it can be taken to be an involution so 
that, by Fact [l], the reversing symmetry group TZ(M) C GL(n, Z) is a semi- 
direct product. Note that if M e GL(n, M) is reversible, it has been shown 
that there always exists an involutory reversing symmetry [|36| , Thm. 2.1]. 
Already for GL(2, Z), this is no longer true [34|: the matrix M = ( ® ^ j is 
reversible in GL(2,Z) with the reversing symmetry G — ( ° t J ) of order 4, 
but with no involutory reversing symmetry in GL(2,Z). 

While reversibility was still a frequent phenomenon in GL(2, Z), see |TIJ, it 
becomes increasingly restrictive with growing n. To see this, let us consider 
necessary conditions for reversibility. If M e GL(n, M.) is reversible, then 
M _1 = GMG^ 1 for some matrix G, and M and M _1 must have the same 
characteristic polynomial, P(x). On the other hand, the spectrum of M 
must then be self-reciprocal, i.e. with A also 1/A must be an eigenvalue, with 
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matching multiplicities. Recall that P(x) = YYi=i( x ~ Aj) an d observe that 

nK)-£SjjG-*)- w 

But by assumption, nr=i( x ~ Aj) = FIILi^ — X")> so we arrive at 

Proposition 2 A necessary condition for the matrix M E GL(n, R) to be 
reversible is spec(M) = spec(M _1 ) and thus the equation 

f _1 \n r n 

™ = kk p(1/l) (23) 

which we call the self-reciprocity of P(x). □ 
One immediate consequence is the following. 

Corollary 4 If M E GL(n, R) «s reversible, det(M) = ±1. J/ ; m addition, 
the multiplicity of the eigenvalue A = — 1 is even (allowing for multiplicity 
i/— 1 is not an eigenvalue of M), then det(M) = 1. 

PROOF: Observe that the reversibility of M implies G = MGM. Taking 
determinants gives the first assertion because neither M nor G is singular. 

Next, note that A = ±1 are the only complex numbers with A = 1/A. 
All other eigenvalues come in reciprocal pairs. Since the determinant is the 
product over all eigenvalues, the second statement follows. □ 

It might be instructive to reformulate Proposition |2] and Corollary f| in 
terms of elementary symmetric polynomials. Let Sk, k — 0, 1, ... , n, denote 
the k-th elementary symmetric polynomial in n indeterminates. The Sk are 
given by So = 1 and 



S k {x 1} ...,X n ) = X h'---' X i k - ( 24 ) 

They are algebraically independent over Z and have the generating function 

n n 

^S h ( Xl ,...,x n )t k = Hil + xJ) (25) 

k=0 i=l 

where t is another indeterminate. 
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Now observe that, for a characteristic polynomial P(x), we have 

n 

P(x) = X n + J2(- 1 ) k Sk(\l,...,\n)x n - k . 

k=l 

Consequently, 

n 

x n P(l/x) = l+^2{-l) k 3 k {X 1 ,...,X n )x k 

k=l 

and a comparison with Proposition ^ reveals that 

S k (X u ... , A n ) = det(M) • S n ^ k (X 1: ... , A n ) (26) 

for all < k < n. For k = 0, this is just the statement that det(M) = ±1. 
Note that the elementary symmetric polynomials, when evaluated at the 
roots of P(x), reproduce (up to a sign) the entries of the last row of the left 
companion matrix (|5p. 

Turning now to the reversibility of matrices in GL(n, Z), we first observe 
that, generically, the reversible cases can only occur when n is even and 
det(M) = +1: 

Proposition 3 Consider M G GL(n, Z) and let P(x) be the characteristic 
polynomial of M. If n > 1 is odd or det(M) = —1 we have: 

(a) if M is reversible in GL(n,Z) ; P(x) is reducible overZ, and the spec- 
trum of M contains 1 or —1; 

(b) if P(x) is irreducible overZ, M cannot be reversible in GL(n, Z). 

PROOF: If M is reversible, A e spec(M) implies 1/A G spec(M), so the 
eigenvalues are either ±1 or have to come in pairs, A ^ 1/A. If n is odd, we 
must have at least one eigenvalue that is ±1, and that gives a factor (x =F 1) 
in P(x). On the other hand, if det(M) = —1, we must have at least one 
eigenvalue that is —1 which gives a factor (x + 1) in P(x). In both cases, 
one notes that whenever ±1 is a zero of a polynomial over Z, factoring out 
(x =F 1) can be done over Z. 

Conversely, if P(x) is irreducible over Z, spec(M) cannot contain an eigen- 
value of the form ±1, and n odd or det(M) = — 1 is then incompatible with 
M being reversible. □ 
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It is clear from this that reversibility is rather restrictive. If P(x) splits 
into irreducible components Pi(x), then each is subject to the constraints 
described above, or has to be matched with its reciprocal partner polynomial 
- if that would be an integer polynomial at all. In particular, if P(x) is 
reducible but contains an isolated irreducible factor of odd order > 3, or 
of even order with constant term —1, reversibility of M is ruled out. For 
example, this confirms that Mi and M2 of ([15]) are not reversible (in fact, 
they are not even reversible in GL(3, R)). 

The key problem in deciding upon similarity of M and M~ l in GL(n, Z) is 
that Z is not a field. But it is clear that the corresponding similarity within 
GL(n, Q) (the matrix entries now belonging to the field of rationals) is both a 
necessary condition and a much easier problem. It would not help to further 



extend Q to R due to the following result, see p3| , Cor. XIV. 2. 3]. 

Fact 8 A matrix M e GL(n, Z) is similar to M~ l within the group GL(n, R) 
if and only if this is already the case in GL(n, Q). □ 

In the light of this, let us first recall some facts about normal forms over 
Q, where similarity is (in theory) a decideable problem. The normal form 
of a matrix M is based on its polynomial invariants, or invariants for short, 
see [[23|, Sec. XIV. 2]. They are often also called the invariant factors of M 
(or, more explicitly, of (xl — M)), compare [^|, Def. 4.4.6], meaning certain 
polynomials that derive from the matrix (xl — M), see below. The following 
result is a direct consequence of Thm. XIV.2.6] or |3], Thm. 5.3.3]. 

Fact 9 Two matrices in Mat(n, Q) are similar in GL(n, Q) if and only if 
they have the same polynomial invariants. In particular, this applies to M 
and M~ l for any M G GL(n, Z) . □ 

Let us briefly recall how the polynomial invariants q l ,...,q r of a matrix 
M e Mat(n, Z) can be found, where r < n is a uniquely determined integer 
that depends on M. We formulate this for integer matrices, but it applies, 
with little change, also to rational ones. Set p = 1 and let p k (for 1 < k < n) 
be the greatest common divisor of all minors of (xl — M) of order k, so 
that p k clearly divides and p n = P(x) = det(xl — M). Let £ denote 

the largest integer k for which p k = 1 and define q i = p e+i /p i+i _!, where 
\<i<r = n — £. These polynomials over Z are the polynomial invariants 
of M and satisfy the following divisibility property: 

1i I Qi+i ■ (27) 
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The prime factors of q i over Z, taken with their multiplicity, are called its 
elementary divisors. The product of the invariant factors of M (equivalently, 
the product of all their elementary divisors) gives the characteristic polyno- 
mial P(x) of M. Furthermore, the minimum polynomial Q(x) of M is given 
by q r , or, equivalently, by the characteristic polynomial P(x) divided by p n _ l . 

Note that a systematic way to find the invariant factors of M is to bring 
(xl — M), seen as a matrix over the principal ideal domain Z[x], into its 
so-called Smith normal form, see || Ch. 5.3] for details. This is a diagonal 
matrix in of the form diag(l, ... , 1, q^x), ... , q r (x)). For large n, calculating 
this form can be a computationally difficult exercise; for small n, the Smith 
normal form can be found from algebraic program packages. Nevertheless, 
significantly, the invariant factors completely determine the Frobenius normal 
form of the matrix M: 

Fact 10 Let M G Mat(n, Z) have polynomial invariants q 1 , ... , q r of degrees 
n 1 , ... , n r , with n l + . . . +n r — n. Then M is similar, in GL(n, Q) ; to a block 
diagonal matrix [E>i, ... , B r ] where E>i is the x rij left companion matrix of 
the polynomial q { . □ 

The existence of a block diagonal matrix similar to M is equivalent to 
the statement that M leaves invariant a set of (cyclic) subspaces of Q™ with 
respective dimensions [ P3"| , Thm. XIV. 2.1] for details. One can 

actually give more refined normal forms by using the elementary divisors of 
each invariant to replace the diagonal blocks Bi with subblock decompositions 
based upon the elementary divisors and their multiplicities. Combining Fact 
|9] and Fact |1(], matrices with the same polynomial invariants can both be 
brought to the same normal form and thus are similar. 



The normal form of Fact [10] highlights the left companion matrices Bi. 
For what follows, we are interested in the reversibility of such matrices. In 
this respect, let and be the left and right companion matrices cor- 
responding to a polynomial P(x). Suppose P(x) conforms to the reciprocity 
condition (23) of Proposition ||. Then one can check that is the in- 



verse of M( £ \ But we already know from (||) that = RM^R~ l where 
R = R~ l is the involution from (0). Combining this with the normal form 
above, we obtain: 

Theorem 2 Let M 6 GL(n, Z). Then, M is reversible in GL(n, Q) if and 
only if each of the polynomial invariants of M satisfies the reciprocity condi- 
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tion (|23| ) separately. In this situation, the reversing symmetry can be chosen 
to be an involution. 



Proof: By Fact [TO, M = SDS' 1 where S G GL(n, Q) and D G GL(n, Z) 



is a block diagonal matrix of the form D = [B\, ... ,B r ], where r > 1 and 
-Bj G GL(n.j, Z) is the left companion matrix corresponding to the invariant 
q { of degree m. It follows that M' 1 = SD^S' 1 where D' 1 = [B^ 1 , ... , B' 1 ]. 
Consequently, M and M _1 are similar if and only if D and D~ l are similar. 

Suppose that each of the polynomial invariants satisfies the condition (p3[) . 
Then, from the remark before Theorem [| £?~ is the right companion matrix 
corresponding to q { and is similar to Bi via the involution Ri G GL(nj,Z) 
which consists of l's on its ant i- diagonal, as in (|10|). It follows that D is 
similar to D~ l via the block diagonal involution R := [R 1 , ... , R r ], and so M 
and M _1 are similar by the involution SRS^ 1 . 

On the other hand, suppose that M is similar to M~ l in GL(n, Q). Hence, 
the corresponding block diagonal matrices D and D~ l are also similar in 
GL(n,Q), via some element G. Now, to each block Bi of D corresponds 
an invariant vector subspace Vi of Q n of dimension rij. A subspace is 
thus either mapped by G to itself (it is a symmetric subspace) or to another 
subspace Vj of the same dimension. In the first case, Bi must be conjugate to 
its inverse via the restriction of G to Vi. This means that the characteristic 
polynomial of Bi, which is q i: must satisfy the condition On the other 
hand, if Vi is mapped to Vj by G with rij = nj, it follows that the invariants 
q i and g. differ by at most a sign. They thus share the same eigenvalues and 
must each satisfy condition ( p3|) on their own. □ 

If M has only one non-trivial invariant, g 1 (x), it follows from the above 
discussion that its characteristic polynomial P(x) coincides with its minimal 
polynomial Q(x) and both equal q t (x) (so M is cyclic from Fact §). Con- 
versely, M cyclic means it has only one invariant. The previous Theorem 
now gives: 

Corollary 5 If M G GL(n, Z) has only one polynomial invariant, in partic- 
ular if the characteristic polynomial P(x) is irreducible over Z, then M is 
reversible in GL(n, Q) if and only if its characteristic polynomial P{x) sat- 
isfies the reciprocity condition (E3). □ 



Note that Theorem ^| and Corollary [5] do not extend to requiring, for 
reversible M, that the elementary divisors within an invariant polynomial 
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satisfy fl23|). For example, any matrix in GL(4, Z) with one invariant poly- 
nomial qi(x) = P(x) = (x 2 — x — l)(x 2 + x — 1) is reversible in GL(4, Q) 
although the elementary divisors separately violate (^). 

Let us give some illustrations of the use of Theorem ^| and Corollary |5] for 
small values of n. These results show that all M G SL(2, Z) are reversible in 
GL(2, Q) because their invariant factors fall into one of the following cases: 

1. qi(x) = q 2 (x) — (x ± 1), so r = 2 and M = =pl; 

2. qi(x) = P(x) = x 2 — tr(M)x + 1, so r = 1 and P(x) is self-reciprocal. 



Yet, we know from [HJ that SL(2,Z) matrices exist that are not reversible 
in GL(2,Z) (see also Section |0| below for further discussion). Furthermore, 
if M G GL(2,Z) with det M = —1, then it can only have one polynomial 
invariant, qi{x) = P(x) = x 2 — tr(M)x — 1. By Proposition § or Corollary ||, 
M is reversible in GL(n, Z) if and only if M has eigenvalues A = ±1 and qi(x) 
factors into (x — l)(x + 1), in which case M G GL(2, Z) is an involution and 
reversible, with reversing symmetry as itself. This approach gives another 
way of retrieving some of the results of [|1(J on the reversibility in GL(2,Z). 

Turning to GL(3,Z), Proposition ^| implies that if M is reversible then 
P(x) must have a factor (xil). In other words, (x±l) 1 , for some 1 < i < 3, is 
an elementary divisor of P(x). If M has more than one polynomial invariant, 
then the divisibility property ( |2"7| ) implies that P(x) completely decomposes 
into a product of 3 factors of the form (x± 1). Generically, however, this will 
not happen and instead M is reversible in GL(3, Q) if and only if it has only 
one invariant of the form q x (x) = P(x) = (x ± l)(x 2 — (tr(M) ± l)x + 1). 

For GL(4,Z), elements with three or four polynomial invariants are di- 
agonal matrices with +l's or (an even number of) — l's on the diagonal. 
They are all reversible. Reversible elements with two invariants must have 
Qii x ) = (x±l) and q 2 (x) = (x ± l)(x 2 — (tr(M) ± 2)x + 1), or q^x) = q 2 {x), 
a monic quadratic with constant term +1. 

As n increases, Theorem ^| can exclude many matrices from being re- 
versible in GL(n, Z) because they are not reversible in GL(n, Q). In particu- 
lar, we can ask for an example M with more than one polynomial invariant 
where the characteristic polynomial P(x) satisfies (p3l), yet M is irreversible 
in GL(n, Q) because it violates Theorem [| If we take n > 2 and even, the 
first possibility appears in GL(8, Z). For example, we can take a matrix with 
invariants q^{x) = x 2 — x — 1 and q 2 {x) = (x 2 — x — l)(x 2 + x — l) 2 . Both 
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polynomials violate the reciprocity condition (p3| ) but they compensate each 
other so that their product, the characteristic polynomial, does satisfy the 
condition. The Frobenius normal form with these invariants is the block 
diagonal matrix M = [Bx, B 2 ], 



M 



/o 

1 






v° 



1 
1 























1 






-1 







1 




-4 








1 




3 









1 


4 



o\ 







1 

-v 



(28) 



This matrix is not reversible in GL(8,Z) (although its square is, see Section 
U below). 

The remaining problem is now to find possible reversibility in GL(n, Z) 
of unimodular matrices that are already reversible in GL(n, Q). For n = 2, 
we were able to solve this problem using a special algebraic structure (the 
amalgamated free product) of PGL(2, Z). This structure is not available for 
n > 3. 

Important examples of unimodular integer matrices which are reversible 
in GL(n, Q) are the symplectic matrices Sp(2n, Z) C SL(2n, Z). Recall that 
a symplectic matrix M G Sp(2n, M) satisfies M l JM = J where M l denotes 
the transpose of M and J is the 2n x 2n integer block matrix 



J 



of order 4. Since in general 

M'.IM = J 



1 
-1 



JM~ l J~ 



(29) 



it follows that M G Sp(2n,lR) is reversible if and only if M is similar to 
M* in GL(2n,IR). But any invertible square matrix with entries in a field 
F is similar to its transpose in GL(m,F), see ]3], Prop. 5.3.7] (but this need 
not be true e.g. in GL(m, Z)). In particular, M G Sp(2n, Z) is reversible in 
GL(2n, Q) and its invariant factors will all satisfy (p3|). Also, it is clear from 
( |2"P| ) that if M is symplectic and symmetric, then M l = M = JM^J^ 1 so 
that M is actually reversible in GL(2n, Z). 
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Orthogonal integer matrices U G GL(n, Z) which satisfy [/[/* = 1 are 
other examples of unimodular integer matrices which are always reversible in 
GL(n, Q) (since U^ 1 = U l and, as before, U* and U are similar in GL(n, Q)). 

We make some further remarks on the problem of deciding reversibility 
within GL(n, Z). Firstly, note that in the proof of Theorem Q above, when 
we used the reversibility of the left companion matrices with characteristic 
polynomials satisfying (0), this reversibility was in GL(n, Z) itself. Further- 
more, the reversing symmetry R was an involution, so by Fact [l| we have the 
following result. 

Theorem 3 For each integer polynomial P(x) of degree n that satisfies the 
necessary self-reciprocity condition ( |2"3"D for reversibility, there is at least one 
reversible matrix class in GL(n, Z), represented by the left companion matrix 
MW ofP(x), and we have TZ(M^) = S(M®) x s C 2 . □ 

Secondly, by Fact |3|, the number of representing matrix classes of an irre- 
ducible characteristic polynomial P(x) equals the class number of the order 
Z[a], with a any of the roots of P(x). If this class number is one, there is 
only the class represented by the companion matrix. If the class number is 
two, one is the companion matrix class which we know to be reversible if the 
spectrum is self-reciprocal. But then, the other class must also be reversible 
because there is no further partner left. 

Class numbers are widely studied in algebraic number theory, and one 
can find both extensive tables in books (e.g. see |17| , |30f ) and also various 
program packages to calculate them, e.g. the program package KANTQ. 

Let us add another example, of rather different flavour, and look at an 
interesting class of algebraic integers, the so-called Salem numbers. They are 
the algebraic integers a > 1 with all conjugates a' having modulus \a'\ < 1 
and with at least one conjugate on the unit circle, see |35], Ch. III. 3] for 
details. Salem's Theorem then says that their degree is always even, that a 
and 1/a are the only real conjugates, and that all other conjugates are on 
the unit circle. In particular, a Salem number is a unit. Putting our above 
results to work, we get 

Corollary 6 Each Salem number occurs as the eigenvalue of a reversible 
toral automorphism. If a is a Salem number of degree n = 2m, and M a 
corresponding GL(n, Z) matrix, then S(M) ~ C 2 x Z m . □ 

7 Sec |http: //www, math. TU-Berlin.de/^kant/| . 
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The polynomial P(x) = x 4 — 2x 3 — 2x 2 — 2x + 1 provides one of the 
simplest examples. Its roots are r ± yjr (both real) and — (r — 1) ± i\Jr — 1 
(both on the unit circle), where r = (l + v^o") /2 is the golden number. The 
corresponding companion matrix M is not of finite order, and the reversing 
symmetry group is thus TZ(M) ~ (C 2 x I?) x s C 2 . 

Let us come back to the general discussion and ask for the properties 
of reversing symmetries. Let G be a reversing symmetry of M, so that 
GM = M~ l G and hence also GM n = M~ n G for all n G Z. If p(x) is any 
polynomial, we then also get Gp(M) = p(M~ 1 )G. Since G is a reversing 
symmetry, G 2 is a symmetry. If we now assume that M G GL(n, Z) has 
minimal polynomial of degree n, we get G 2 G Q[M] from Fact i.e. G 2 = 
g(M) for some q with coefficients in Q. Consequently, Gq(M)G~ 1 = g(M _1 ), 
but also Gq(M)G~ 1 = G 2 = q(M), so that g(M) = q(M~ 1 ). 

If g(M) is a monomial, i.e. q(M) = M e for some £, then M 2 ^ = 1 and hence 
G A = 1. This case is also discussed in |L6], Prop. 2(i)]. It clearly extends to 
the situation that g(M _1 ) = (q(M)) , which is more general. Apart from 



this, we recall the following result from ]T6|, p. 21] (its proof, which was only 



contained in the preprint version of |16|], is a coset counting argument). 



Fact 11 Let M be of infinite order. If the factor group S(M) / (M) is finite, 
then any reversing symmetry G of M must be of finite order, and G 2k = 1 
for some integer k that divides the order of the factor group. □ 

Let us only add that, in line with our above argument, one first obtains 
G 2h G (M) and hence G 4fc = 1. But (M) ~ Z by assumption, so it cannot 
have a subgroup of order 2, and thus already G 2k = 1. 

Fact [11] certainly applies to our scenario whenever M is not of finite order, 
but S(M) Abelian and or rank 1. This type of result is helpful because it 
restricts the search for reversing symmetries to one among elements of finite 
order. It is certainly possible to extend the result to other cases, but in 
general one has to expect reversing symmetries of infinite order, in particular 
if the rank of S(M) is > 2. Even then some results are possible because it 
would be sufficient to know whether reversibility implied the existence of 
some reversing symmetries of finite order. However, this question is more 
involved and thus postponed. 
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3 Extensions and further directions 



In this Section, we will summarize some additional aspects of our analy- 
sis, namely the extension of symmetries to affine mappings, the modifica- 
tions needed to treat the related situation of the projective matrix group 
PGL(n, Z), and the extension of (reversing) symmetries from a group setting 
to that of (matrix) rings or semi-groups. 



3.1 Extension to affine transformations 

So far, we have mainly discussed linear transformations (w.r.t. the torus), 
but it is an interesting question what happens if one extends the search for 
(reversing) symmetries to the group of affine transformations. Since both 
arguments and results are the exact analogues of those for the case n = 2 as 



derived in [fL0| , we will be very brief here. 

In Euclidean n-space, the group of affine transformations is the semi-direct 
product Q a = R n x s GL(n, R), with R n being the normal subgroup. Elements 
are written as (t, M) with t G R n and M G GL(n, R), and the product of two 
transformations is (t, M) ■ (f, M') = (t + Mt', MM'). The neutral element is 
(0, 1 ), and we have (t, M)" 1 = (-M~H, M' 1 ). 

If we now observe that T n = R n /Z n , it is immediately clear that the affine 
transformations of T n form the group 

gj n = T n x s GL{n,Z) (30) 

which is still a semi-direct product. Here, T n can be written as [0, 1)™ with 
addition mod. 1, and the product of transformations is modified accordingly. 

If we now ask for an affine (reversing) symmetry of a matrix M (now being 
identified with the element (0, M) G Gj") we find 

Proposition 4 The affine transformation (t, G) is a (reversing) symmetry 
of the toral automorphism (0, M) if and only if 

(a) G is a (reversing) symmetry of M in GL(n, Z) and 

(b) Mt = t (mod 1). 

Proof: We have (*, G) ■ (0, M) = (t, GM) and also (0, M ±v ) ■ (t, G) = 
(M^t, M ±X G). But then, the statement follows from the uniqueness of fac- 
torization in semi-direct products. □ 
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From the condition Mt = t (mod 1) it is clear that we need not consider 
all translations in T n but only those with rational components, which we 
denote as A^. For many concrete problems, it would actually be even more 
appropriate to restrict to discrete sublattices, e.g. to the so-called g-division 
points A q ~ (Cg) n which consists of all rational points with denominator q. 
We will not follow this idea here, however. 

From the above result, it is clear that we can get (reversing) fc-symmetries 
(recall the definitions from Section O). In fact, the equation M k t = t on 



the torus has a k = | det(M fe — 1)| different solutions provided no eigenvalue 
of M k is 1. Clearly, 

a k = Y. i ' c i ( 31 ) 

£\k 

where c e counts the true orbits of length £, and the Mobius inversion formula 
gives 

c k = tit) -at (32) 

£\k 

with the Mobius function /x(m) |H| p. 29]. If c k is positive for some k, we 
get a fc-symmetry (and, hence, eventually a reversing /c-symmetry) of M. 
These numbers can easily be calculated explicitly, where a very natural tool 



is provided by the so-called dynamical or Artin-Mazur £- functions |T5|] . Here, 
the a fc 's can be extracted from the series expansion of the logarithm of the 
^-function, while the c fc 's appear as exponents of the factors of the Euler 
product expansion of the ^-function itself. 



3.2 The case of PGL(n, Z) 

Let us start by the observation that PGL(n, Z) can be described via quoti- 
enting w.r.t. {±1}, i.e. 

PGL(n, Z) ~ GL(n, Z) /{± 1} . 

In other words, rather than consider single matrices M, one has to consider 
pairs, [M] := {±M}. Let us write S[M] for the new PGL case and keep the 
old notation for the GL situation treated above. 

The modification needed for the symmetry analysis given above is then 
actually fairly trivial, as we always had ±1 among them, and we can simply 
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factor that out. So, we get 



S[M] ~ S(M)/{±1}. (33) 

The case of reversing symmetries, however, requires some care. Since we 
now calculate mod ±1, a projective matrix [M\ can also be reversible through 
GMG- 1 = -M -1 . But if this happens, the square, M , is again reversible 
in the old sense. This mechanism can (and will) give rise to reversing 2- 
symmetries (recall Section |1.2| for the definition) in GL(n, Z). Whereas [T(J 



gave examples in GL(2,Z), Eq. (^) shows that skew-symmetric symplectic 
matrices satisfying M* = — M are not reversible in GL(2n,Z) whereas their 
squares are reversible. Also, the example (^) is reversible in PGL(8, Z) since 
a G G GL(8, Z) can be found, by direct calculation, that satisfies the relation 
GMG- 1 = -M- 1 . 

Let us finally check what happens in the extension to affine transforma- 
tions. In complete analogy to the case n = 2, see ||10|| , one can show that 
the corresponding affine group is the semidirect product A 2 x s PGL(n,Z) 
with A2 the 2-division points. This really is the consequence of identifying x 
with —x on T n , and A 2 is the set of 2 n translations that satisfy the condition 
t = -t (mod 1). 

Now, the above Proposition (| applies to the case of PGL-matrices in very 
much the same way, just the possible translations t are restricted to the 
2-division points. 



3.3 Symmetries among general integer matrices 

For most of this article, we have focused on matrices in GL(n, Z) and their 
symmetries within the same group. However, none of the proofs given above 
depends on that restriction, and one can indeed also treat the case that both 
M and its symmetries are allowed to live in the larger set Mat(n, Z) which 
is no longer a group w.r.t. multiplication, but a ring. Nevertheless, we will 
continue to use the symbol S(M), now meaning 

S(M) := {G G Mat(n, Z) | [M, G] = 0} . 

The most obvious extended symmetries which one gets in Mat(n, Z) are the 
integer multiples of the identity, but there really is a hierarchy of objects to 
look at, and it is most transparent if one phrases the situation for a matrix 
in Mat(n,Q) first: 
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Fact 12 Let M G Mat(n, Q) and let its characteristic polynomial P(x) be 
irreducible overQ. Let A be any of its roots and K = Q(A) the corresponding 
algebraic number field. Then the following statements hold. 



(a) 


Cen ^Mat(n,« 




~ K. 




(b) 


Cen ^GL(n,Q 




~ K* , where K* = K \ {0}. 




(c) 


CeI1 ^Mat(n,; 


i){M) 


~ O , where O is an order in K . 




(d) 


CeI1 ^GL(n,Z 


)(M) 




□ 



The proof is a slight variation of what we did for Fact [| and Lemma [1|, and 
need not be spelled out again. It is important to note that the order O 
appearing here, as mentioned before, in general is not the maximal order of 
K, though it contains Z[A]. The following is now an immediate consequence. 



Proposition 5 Let M be an integer matrix with irreducible characteristic 
polynomial P(x). Let X be a root of P(x), and let O max be the maximal order 
in Q(A). Then, S(M) is both a Z-module and a ring, and isomorphic to an 
order O that satisfies Z[A] C0C C max . □ 

Note that there is now also a natural extension to the case of simple matrices 
M, compare Theorem [I] and Eq. (0), but we omit further details here. Also, 
from Fact |5| it is clear that S(M) is Abelian if and only if the minimal 
polynomial of M has degree n. 

As to reversibility, this new point of view requires some thought. By 
Corollary [|, M reversible implies det(M) = ±1, so reversible integer matrices 
are restricted to GL(n, Z). It would then not be unnatural to also insist on 
the existence of at least one unimodular matrix G with M _1 = GMG -1 , and 
reversibility is basically as above, except that, if we enlarge the symmetries of 
M from subgroups of GL(n, Z) to subrings of Mat(n, Z) , reversing symmetries 
get enlarged accordingly. Note, however, that the ring structure is lost: the 
sum of a symmetry and a reversing symmetry is not a meaningful operation 
in this context. Together, they only form a monoid, i.e. a semi-group with 
unit element. 

To go one step further, one could then also rewrite the reversibility con- 
dition as 

G = MGM 
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and only demand that G is non-singular, to avoid pathologies with projec- 
tions to subspaces and to keep the statement of Corollary Note that this 
is a slightly weaker form of reversibility, as one does not assume that G^ 1 
is a meaningful mapping in this context. In particular, it is clear that there 
is then no need any more to restrict G to unimodular integer matrices, so 
that now G G GL(n, Q). This will, in general, lead to new cases of (weak) 



reversibility, as is apparent from the explicit constructions in |T0 
To give a concrete example, consider the matrices 

*-(}•) - -g_y. (34, 

Then M is the matrix from [1^, Ex. 2] that was shown to be irreversible in 
GL(2, Z). In fact, the automorphism of T 2 induced by M is irreversible even 
in the larger group of homeomorphisms of the 2-torus, see and |T|, p. 9] 
for details on the connection between general and linear homeomorphisms. 
Nevertheless, one can check that G = MGM, where det(G) = —9. In other 
words, M is reversible in GL(2,Q), as are all elements of SL(2, Z) by our 
previous discussion following Corollary ||. 

Note that G does not induce a homeomorphism of the 2-torus because 
G^ 1 is not an integer matrix. However, G does induce an automorphism on 
any lattice of the torus of the form 

A, := {(? zy\0<m,n<q} 



for which det(G) (mod q). This is relevant as a recent study [|T^] shows: 
the quantum map which corresponds to M (which, in turn, corresponds to the 
action of M on a (Wigner) lattice of the torus) showed an eigenvalue statistics 
according to the circular orthogonal ensemble (COE) rather than the unitary 
one (CUE). So, even though M does not have a reversing symmetry in the 



sense of Section fO| the presence of "pseudo-symmetries" such as G still 
leave their mark! 
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